Implementing topological quantum manipulation with superconducting circuits 
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A two-component fermion model with conventional two-body interactions was recently shown to have any- 
onic excitations. We here propose a scheme to physically implement this model by transforming each chain 
of two two-component fermions to the two capacitively coupled chains of superconducting devices. In par- 
ticular, we elaborate how to achieve the wanted operations to create and manipulate the topological quantum 
states, providing an experimentally feasible scenario to access the topological memory and to build the anyonic 
interferometry. 
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Topological ordered states emerge as a new kind of states 
of quantum matter beyond the description of conventional 
Landau's theory [1], whose excitations are anyons satisfy- 
ing fractional statistics. A paradigmatic system for the exis- 
tence of anyons is a kind of so-called fractional quantum Hall 
states [2]. Alternatively, artificial spin lattice models are also 
promising for observing these exotic excitations iQlUlBl]. 
taev models IHQ] are most famous for demonstrating anyonic 
interferometry and braiding operations for topological quan- 
tum computation. 

Since anyons have not been directly observed experimen- 
tally, a focus at present is to experimentally demonstrate the 
topological nature of these states. Abelian anyons maybe 
relatively easy to achieve and to manipulate in comparison 
with the nonabelian ones, thus it is of current interest to ex- 
plore them both theoretically and experimentally. Kitaev con- 
structed an artificial spin model [3], i.e., the toric code model, 
which supports the abelian anyon. But the wanted four-body 
interactions are notoriously hard to generate experimentally 
in a controllable fashion. Alternatively, it was proposed yfl to 
generate dynamically the ground state and the excitations of 
the model Hamiltonian, instead of direct ground-state cool- 
ing, to simulate the anyonic interferometry. On the other 
hand, implementation of another Kitaev's honeycomb model 
J3] was also suggested in the context of ultracold atoms |6], 
polar molecules [7], and superconducting circuits JUl. The 
honeycomb model H] is an anisotropic spin model with three 
types of nearest-neighbor two-body interactions, which sup- 
port both abelian and nonabelian anyons. It was shown [4] 
that the toric code model can be obtained from the limiting 
case of the honeycomb model. Using this map, preliminary 
operations for topological quantum memory and computation 
were also addressed StlSllIltlll]- Nevertheless, in this case, 
anyons are created by the fourth-order perturbation treatment, 
which would blur the extracted anyonic information |[l2tl . In 
addition, this map is good but not exact, so that one may get 
both anyonic and fermionic excitations. Therefore, new meth- 
ods for implementing the model and manipulating the relevant 
topological states are still desirably awaited. 



Recently, a two-component fermion model [13] with con- 
ventional two-body interactions was shown to have anyonic 
excitations, which obey the same fusion rules as those of the 
toric code model and are mutual semions. This model is 
promising because it provides an example for abelian anyons, 
which can be directly implemented. In this Rapid Communi- 
cation, we propose to physically implement this model with 
appropriately designed superconducting circuits. In particu- 
lar, we elaborate how to achieve the wanted operations that 
create and manipulate the topological states as well as anyons 
with the cavity-assisted interactions using an external mag- 
netic drive, providing an experimentally feasible scenario to 
access the topological memory and to build the anyonic inter- 
ferometry. 

The Hamiltonian for the two-component fermion model in 
a two-dimensional square lattice, as shown in Fig. [Tf a X is lfl3ll 

Hfi = -J (? ^(2n T , i -l)(2n TJ -l) 

-J p ^2(2n Ui -l)(2n u -l) 

(id) 

+f/^(2n T , i -l)(2n;, i -l), (1) 



where n. 



are annihilation operators of spin-s 
fermions, and (i, j) mean the nearest neighbors along the hor- 
izontal diagonals of squares. The ground states of this Hamil- 
tonian are highly degenerated, i.e., every individual chain is 
ferromagnetic. As shown in Ref. Il3ll . the low-lying excita- 
tions are deconfined mutual semions under the open boundary 
condition. 

To show the nature of the low-lying excitations, we map 
the square lattice to the honeycomb lattice, as shown in Fig. 
[TJb), by extending each lattice site in the square lattice to be a 
vertical link of Majorana fermions [see Fig. [He)] defined by 

i% = -*( c T,i- c T,i)'^ = c T,i + c |,i ; ^ = -*( c i^~ c L)' 
and Xiu, = c l,i + c \ i- These "real" fermion operators obey 
= Xi =1- Otherwise, they are anticommutative. In 
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FIG. 1 : (Color online) A map from a square lattice to a honeycomb 
lattice by extending each square lattice site to be a vertical link, (a) 
The square lattice, where denotes the nearest neighbors along 
the horizontal diagonal direction, (b) The honeycomb lattice, where 
ib, w label the black and white sublattices. The indicated zig-zag 
chain is one of the lines for the Jordan- Wigner transformation, (c) 
Each square lattice site being extended as a vertical link of Majorana 
fermions with i w (it) labels the white (black) sublattice and i w on 
the top of ii,. (d) Site-labels within a plaquette. 



this Majorana fermion representation, the Hamiltonian (Q]) is 
mapped to 

(*i> (ij) i 



where W, 



iQj with Qi = iijji w il>i b and W l . 



with Qi = ixi m Xh ■ Since there is no coupling between the 
chains, the present model may be transferred to a spin model 



#/3 = -J q E Wp - J p E - U H s t s t 



(3) 



QZ 

>-3'.<h a 3V 



by using Jordan- Wigner transformation 

to. = s f IL;. , % *h = SJ b U, 
X jw = SlU fs<jw Sp t , and Xjb = SIR, 
where S x ' v,z are the corresponding Pauli matrices for the 
defined Majorana fermions, W P = SfS$S§SlS%SI and 
W P = SfS%S§S%S%SI with the site labels within a plaque- 
tte depicted in Fig. \V[d). The order of the sites is defined 
as follows: j s > j t if the zigzag line [one of such lines is 
indicated in Fig. 02b)] including j s is higher than that of jt, 
or if j s is on the right hand of j t when they are in the same 
line. It is straightforward to check [Wp, Wp>] = for every 
plaquette, and all of them also commute with the third term. 
In this spin model, the ground state can be written as 



\G) = l[(l + Wp)(l + Wp)\ct>), 



(4) 



where \<f>) = |1 • • • 1) 
the eigenvalue of Sf 



is a reference state and each "1" means 
( being 1. Similar to the Kitaev's hon- 
eycomb model all excitations here may be labelled by two 
quantum numbers Wp and Wp. The low energy excitations 
fall into two closed subsets, each can be graded by a Z2 x Z2 
group. The fusion rules of these excitations are equivalent 
to the excitations in the toric code model OL and different 
graded vortices are mutual abelian semions IU5I1 . 




FIG. 2: (Color online) A schematic circuit of two chains of capaci- 
tively coupled superconducting devices, labeled by a and b, to imple- 
ment a chain of two-component fermions. Here only the first device 
of the two chains is explicitly shown, while others are simply denoted 
as the filled circles with different colors (shades) for different chains. 

We now proceed to implement the model with capacitively 
coupled superconducting devices, i.e., the Cooper pair box. 
The key idea is to use two chains of capacitively coupled su- 
perconducting devices, as shown in Fig. |2j to implement a 
chain of two-component fermions. A building block of our 
implementation, as shown in the rectangle of Fig. |2] is the 
two capacitively coupled superconducting devices. A typical 
design of a Cooper pair box consists of a small superconduct- 
ing island with n excess Cooper pair charges connected by a 
Josephson junction with coupling energy Ej and capacitance 
C,j. A control gate voltage V g is applied via a gate capacitor 
C g . To quantize the circuit equation, we first introduce the 
Hamiltonian and then convert the classical momentum vari- 
able to the momentum operator. Then the Hamiltonian reads 



^m 2 



E 1 } cos ip" 



a , • b 



— adp a (p 



(5) 



where tp* 1 is the gauge phase drop across the corresponding 



junction, C? = C% + C C with C% = C'J + C% a = (h/2e) 
and the induced charge = C^V* 1 /2e. At temperatures 
much lower than the single-pair charging energy, i.e., kpT <C 
E% = e 2 / (2Cq ), and restricting the gate charge to the range 
of n„ G [0, 1], only a pair of adjacent charge states {|0), |1)} 
on the island are relevant. The Hamiltonian © is then reduced 
to [161 



H, 



q2 



5S> 



where e a ^ = 

A" = E 1 }, A = 



1 



6(a) 



2< 



a2 + AV2] + \a a z a b z , (6) 



C c )/A with A = CfC b t - Cl 



2e 2 (C, 

e 2 C c /A, T\ £ {a, 6} and a XiZ denotes the cor- 
responding Pauli matrix in the basis of {|0), |1)}. The single- 
device terms in Hamiltonian © can be tuned to be zero by 
conventional methods B17I1 . Therefore, in what follows, we 
do not take them into consideration. For two identical devices 
(Cg = C ), e = AE C and A = e 2 C c /(C 2 - C 2 ) ~ 2[3E C with 
f3 = C c /Cq. It is notable that the strength of this interaction, 
proportional to the coupling capacitance, is stronger than any 
other present-known coupling methods. 
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The circuit Hamiltonian of the two coupled chains, as de- 
picetd in Fig. [2j can be obtained in a similar way. In this 
extended multipartite coupling case, the long-range interac- 
tion between the devices would appear, which decays expo- 
nentially as Z?' 2 "- 7 ' with i and j being the site labels of the two 
involved devices lfl8ll . Therefore, the long-range interaction 
is negligible if (3 <C 1, and in typical experiments (3 ~ 0.05 
lfl9ll . Up to the first order of (3, the interaction Hamiltonian is 
given by 

hjj = £ w: ;( i + D + £ x ^<v C7) 

V\3 3 

where A, ; ~ e 2 C v /[C + 2(C C + 2C V )} and A c ~ e 2 C c /[C* + 
2(C C + C a + C b )]. 

Once the coupled chains are placed according to the ge- 
ometry of Fig. 0I a X a corresponding two-dimensional square 
lattice model is constructed. Drawing an analogy between the 
device states in chain a (b) and spin f (J.), the above interac- 
tion Hamiltonian in the addressed system may be rewritten as 
the two-component fermion model of Eq. CO, with the pa- 
rameters (A a , Xb, A c ) corresponding to ( J q , J p , U) 0]. As a 
result, the topologically protected ground state of Eq. may 
be implemented with the present setup of superconducting de- 
vices. 

In addition, to accomplish certain topological quantum 
manipulation tasks, including the examination of the anyon 
statistics, it is a must to have a set of basic operations of the 
devices. In the Majorana fermion representation, the wanted 
operations are fl3ll 

Sf b = iXh^h, S z jw = iXj w *Pj w , (8a) 
s l = *h II (iXi1>j), SI = Xju II (iXiti), (8b) 

3<3b j<3w 

where is the spin-flip operator for a given site and it also 
transfers a double fermion occupancy to empty or vice versa. 
SJ denotes a nonlocal operation that creates or annihilates a 
fermion at site j and changes the site occupation for sites j < 
jb(w) Q- In the present model, S z and S x are effective Pauli 
matrices, which, according to Ref. [4], may create and move 
the excitations. 

At this stage, let us elaborate how to obtain the wanted op- 
erations in Eqs. dial ) and d8bl ). Notably, individual address- 
ability is normally a prerequisite in such manipulation. In the 
present proposal, the size of the device setup is macroscopic, 
thus individual addressability is taken as granted. To manip- 
ulate the states, we put the lattice into a microwave cavity, 
with the geometry of the hybrid system being explained in the 
caption of Fig. [3] For simplicity, we consider only the single- 
mode standing wave cavity. To be more specific, we reach the 
cavity assisted manipulation by a magnetic drive l2lll . The 
interaction can be switched on/off by modulating the external 
magnetic field to be ac/dc 1 2 1 ] . With the dc magnetic flux, the 
external flux is merely used to address separately the single- 
qubit rotations. Under the cavity field, it is also readily pos- 
sible to tune off single-qubit terms. The wanted operations in 
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FIG. 3: (Color online) A schematic diagram of the cavity-assisted 
manipulation. The x and y coordinates are denoted by the arrows, 
and the z direction is pointing out the xoy plane. The square lattice 
(red) is placed to be parallel to the yoz plane. The square lattice and 
the auxiliary device (blue rectangle) are placed with the x coordinate 
at the antinodes of the single-mode standing-wave cavity. All super- 
conducting devices are placed with their loop plane being parallel to 
the xoz plane, which is perpendicular to the magnetic component of 
the cavity field, letting it be the only contributed component. 

Eq. d8al ) for selected device can be achieved with the cavity 
mediated integrations by tuning the driven magnetic flux of 
the device to be of ac. The Pauli matrices Sj w and Sj b in Eq. 
d8ab correspond to a x ® a x and a y ® a y two-body interaction 
of two devices of jth site, respectively. These two type inter- 
actions for each lattice site can be directly engineered in our 
hybrid implementation l2lll as it allows selected addressing 
of designated devices. These two interactions are mediated 
by the virtue cavity photon, thus we need to keep the cavity 
mode in the vacuum state. Here, the devices work in their 
degeneracy points. 

The common cavity mode can also be used to realize the 
global stringlike S x operators in Eq. ( fSbt . The off-resonant 
interaction between the cavity mode and the selected devices 
isJH 

Hqnd = xnc^2°], (9) 

3 

where n c = a' a is the photon number operator of the cav- 
ity mode, the coupling strength is \ — 5 2 /2^ with g as the 
single-photon Rabi frequency for the cavity mode, and 8 is 
the detuning between the cavity mode frequency lo c and opti- 
cal transition frequency in atomic spins. In our implementa- 
tion, this can be the |1) — > |2) transition of the selected de- 
vices, where |2) is an ancillary energy level beyond the qubit 
subspace {|0), |1)}, and the frequency of the drive ac flux sat- 
isfies lu = u>i2 +lu c + 8 with = 2E c (3 — 2n g )/h. To avoid 
the transitions |0) — ► |1) and |0) — > |2) by the ac drive, we 
tune cDqi = 2E C (1 — 2n g )/hvia n g so that 5 <C Ai :2 , where 
Ai = lu — ojQi and A2 = cDcn +^12 — w are the corresponding 
detunings. This quantum non-demolition (QND) Hamiltonian 
(0 preserves the photon number of the cavity mode. Within 
the n c € {0,1} subspace, the evolution of the QND Hamilto- 
nian during the interaction time t = u j2\ yields 1 1(111 

f I for n c = 

U = exp[-iHT ] = < N (10) 

{ rij °j f ° r n c = 1 
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where N is the number of the selected devices. From Eq. 
( fTUl i, (controlled) string operations for an arbitrary string can 
be achieved fioll . 

If the cavity is initially prepared in the n c = 1 state, 
the global operation reduces to the string operation U z — 
Y[j (Jj ■ As all string operators are equivalent to U z up to 
local single spin rotations ifloll . all string operations for ar- 
bitrary string can be achieved: U x = Ylj a 'J = HU Z H and 
U y = Ylj a) = RU Z R, where H = Ylj Hj and R = \[. Rj 
with Hj = (<rj + cr| ) /v2 being the Hadamard rotation and 
Rj = exp (— if cr|)- Therefore, with this elementary opera- 
tion, creation and manipulation of anyons are likely feasible 
in our scheme. For example, 

s h = ^ n n ttxk^kj 

kt<jb k w <ji, 

= <j Y[ {i<kA;k b ) II «^<fc„)>(H) 
fed <jb fc lu <jb 

denotes a nonlocal operation that could create a domain-wall- 
like excitation/semion (at the site j) from the ground states 
under the open boundary condition [13]. 

If the cavity is initially prepared in a superposition of zero- 
and one-photon states, the global operation in Eq. < TT~Ob re- 
duces to a controlled-string operation: U cs — /i|0}(0| <E> I + 
f|l)(l| ® U z , where the parameters /i and v are controlled 
by the initially prepared photon number state. With such a 
controlled-string operation, one is able to access the topolog- 
ical memory and to build anyonic interferometry 1 10]. 



In simulating the string operations, we need to engineer 
the cavity number states. Therefore, beside the square lat- 
tice, we also place an ancilla device in the cavity, as shown 
in Fig. [3j which is used to control the cavity photon number 
state by swapping its states with that of the cavity using the 
resonate cavity-device interaction. This swap operation can 
be achieved by the famous Jaynes-Cummings model Hamil- 
tonian: Hjc = f2 {a<j + + aJa), which can be implemented 
in our system by choosing the frequency of the ac driven mag- 
netic flux for the ancillary device satisfying u> = oj c + cjoi- 
In this case, we need to tune the device slightly away from 
the degeneracy point, which results in a shorter decohenrence 
time. Fortunately, the resonate operation is also much faster. 

In summary, we have proposed an exotic scheme to imple- 
ment a two-component fermion model using superconducting 
quantum circuits, which was shown to support abelian any- 
onic excitations. Most intriguingly, we have elaborated how 
to achieve all the wanted operations that could create and ma- 
nipulate the anyonic states. Our approach provides an exper- 
imentally feasible scenario to access the topological memory 
and to build the anyonic interferometry. 

We thank Yue Yu, L. B. Shao and Liang Jiang for 
many helpful discussions. This work was supported by the 
RGC of Hong Kong under Grants No. HKU7045/05P and 
No. HKU7049/07P, the URC fund of HKU, the NSFC 
under Grants No. 10429401, No. 10674049 and No. 
10625416, and the National Basic Research Program of 
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